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Abstract 

The first KamLAND results are in a very good agreement with the predictions 
made on the basis of the solar neutrino data and the LMA realization of the MSW 
mechanism. We perform a combined analysis of the KamLAND (rate, spectrum) and 
the solar neutrino data with a free boron neutrino flux fs- The best fit values of 
neutrino parameters are Am 2 = 7.3 • 10 -5 eV 2 , tan 2 9 = 0.41 and fg = 1.05 with the 
la intervals: Am 2 = (6.2 - 8.4) • 10 -5 eV 2 , tan 2 9 = 0.33 - 0.54. We find the 3d upper 
bounds: Am 2 < 2.8 • 10~ 4 eV 2 and tan 2 9 < 0.84, and the lower bound Am 2 > 4 • 10~ 5 
eV 2 . At 99% C.L. the KamLAND spectral result splits the LMA region into two 
parts with the preferred one at Am 2 < 10 -4 eV 2 . The higher Am 2 region is accepted 
at about 2a level. We show that effects of non-zero 13-mixing, sin 2 #13 < 0.04, are 
small leading to a slight improvement of the fit in higher Am 2 region. In the best 
fit point we predict for SNO: CC/NC = 0.33±g;g| and A S D N N ° = 2.8 ± 0.8% (68% 
C.L.), and A^ff < 9% at the 3<r level. Further improvements in the determination 
of the oscillation parameters are discussed and implications of the solar neutrino and 
KamLAND results are considered. 



1 Introduction 



The first KamLAND results [1] are the last (or almost last) step in resolution of the long- 
standing solar neutrino problem [2]. "In the context of two flavor neutrino oscillations 
with the CPT invariance, the KamLAND results exclude all oscillation solutions of this 
problem but the 'Large mixing angle' solution" [1]. In fact, KamLAND excludes also all 
non-oscillation solutions based on neutrino spin-flip in the magnetic fields of the Sun, on 
the non-standard neutrino interactions, etc.. More precisely, KamLAND excludes them as 
the dominant mechanisms of the solar neutrino conversion. 

Soon after the suggestion of the MSW mechanism [3], the "MSW triangle" [4] had been 
constructed in the Am 2 — sin 2 29 plane which corresponds to 1/4 - 1/3 suppression of the 
Argon production rate [5]. The region around the upper right corner of this triangle, or 
large sin 2 29 part of its horizontal side ( "the adiabatic solution" ) is what we call now the 
LMA region. 

Kamiokande [6], SAGE [7] and GALLEX [8] have disintegrated the triangle into pieces 
excluding some parts of oscillation parameter space. First, the results from Kamiokande 
(the absence of a strong spectrum distortion as well as the absence of strong Day-Night 
effect) have chipped off the LMA region. Then Gallium experiments SAGE and GALLEX, 
confirming LMA, have splitted the diagonal side of the triangle into SMA and LOW. Apart 
from that vacuum oscillation solutions (VO) were always present. 

For a long time SMA was the favored solution and LMA was considered as an non ex- 
cluded possibility. Things changed in 1998 when Super-Kamiokande [9] has testified against 
SMA (the flatness of spectrum and the absence of a peak in the zenith angle distribution 
of events in the Earth core bin). Soon after that, the analysis of the Super-Kamiokande 
data for 708 days of operation allowed us [10] to conclude that solar neutrino results provide 
hints that the LMA solution could be correct (see also analysis in [11]). Since that time 
LMA continuously reinforced its position. 

Super-Kamiokande [9] and then SNO [12] have accomplished genesis of the triangle. 
They have practically excluded SMA, disfavored LOW and VO, and shifted the LMA re- 
gion to larger Am 2 (small Day-Night asymmetry at Super-Kamiokande) and smaller mixings 
(small ratio of CC/NC events at SNO). 

By the time of the KamLAND announcement, the solar neutrino data [5, 7, 8, 9, 12, 
13, 14] have definitely selected LMA as the most favorable solution based on neutrino mass 
and mixing [15, 16]. The best fit point from the free boron neutrino flux fit [16] is 

Am 2 = 6.15 • 10" 5 eV 2 , tan 2 9 = 0.41, f B = 1.05, (1) 

where f B = F B /F§ SM is the boron neutrino flux in the units of the Standard Solar Model 
predicted flux [17]. 
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On basis of the solar neutrino results (and the assumption of the CPT invariance) 
predictions for the KamLAND experiment have been calculated. A significant suppression 
of the signal was expected in the case of the LMA solution. The predicted ratio of the 
numbers of events with the visible (prompt) energies, E p , above 2.6 MeV with and without 
oscillations equals [16]: 

R%Z A = 0.651838 (3a). (2) 

For other solutions of the solar neutrino problem one expected R KL = 0.9 — 1, where the 
deviation from 1 can be due to the effect of nonzero 1-3 mixing. 

In the best fit point (1) the predicted spectrum has (i) a peak at E p « (3.0 — 3.6) MeV, 
(ii) a suppression of the number of events near the threshold energy E p « 2.6 MeV and (iii) 
a significant suppression of the signal (with respect to the no-oscillation case) at the high 
energies: E p > (4 — 5) MeV [16]. No distortion of the spectrum is expected for the other 
solutions. 

The first KamLAND results: both the total number of events and spectrum shape [1], 
are in a very good agreement with predictions: 

R™1 = 0.611 ± 0.094 . (3) 

The spectral data (although not yet precise) reproduce well the features described above. 
As a result, the allowed "island" in the Am 2 — sin 2 26 plane with the best fit KamLAND 
point covers the best fit point from the solar neutrino analysis [1]. 

In this paper we present our analysis of the KamLAND data as well as combined analysis 
of the solar neutrinos and KamLAND results. The impact of KamLAND on the LMA MSW 
solution is studied. We consider implications of the combined (solar neutrinos + KamLAND) 
analysis and make predictions for the future measurements. 

During preparation of this paper, several studies of the first KamLAND results and solar 
neutrino data have been published [18] - [25]. Our conclusions are in a good agreement with 
results of those papers, although there are some differences. Detailed comparison of results 
will be done later. 

2 KamLAND 

In a given energy bin a (a = 1, ....13) the signal at KamLAND is determined by 

r E a +AE r 

N a = Aj2J dE p J dE' p P t F t af(E p , E' p ) , (4) 
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where AE = 0.425 MeV, 



P»=(l-sin 2 2tf 12 sin 2 ^^) (5) 

is the vacuum oscillation survival probability for % reactor situated at the distance from 
KamLAND, is the flux from % reactor, a is the cross-section of vp — > e + n reaction, E p is 
the observed prompt energy, E' p is the true prompt energy, f(E p , E' p ) is the energy resolution 
function, A is the factor which takes into account the fiducial volume, the time of observa- 
tion, etc.. We sum over all reactors contributing appreciable to the flux at KamLAND. 

The suppression factor of the total number of (reactor neutrino) events above certain 
threshold is defined as: 

r, /a 2 2 m iV(Am 2 ,tan 2 6>) 

R K L{Am 2 , tan 2 9) = — ^ ^ , (6) 

where N = J2 a N a , N a is given in (4) and N is the total number of events in the absence 
of oscillations (following KamLAND we will call N and N the rates). 

1). KamLAND spectrum. High threshold. We perform x 2 analysis of the KamLAND 
spectrum, defining 

X% ec = E E ( N a~ N?)a£(N b - N?) , (7) 

a=1...13 6=1. ..13 

where a a b is the covariance matrix in which the systematic uncertainties were propagated 
to the spectral bins. 

We find that for E p > 2.6 MeV the minimum of x 2 pec is achieved for 

Am 2 = 7.34 • 10~ 5 eV 2 , tan 2 9 = 0.453, (8) 

and in this point x 2 '/d.o.f. = 2.92/11. Notice that in contrast with the KamLAND result [1] 
our best fit mixing deviates from the maximal mixing. 

We present in fig. 1 the contours of constant confidence level with respect to the best fit 
point (8) in the (Am 2 — tan 2 9) plane using relation: \ 2 = Xmin + Ax 2 , where A% 2 = 1, 3.84 
and 6.63 for la, 95% and 99% C.L. correspondingly. 

The contours manifest an oscillatory pattern in Am 2 in spite of a strong averaging 
effect which originates from large spread in distances from different reactors. The pattern 
can be described in terms of oscillations with certain effective distance, L e ff, and effective 
oscillation phase e // : 

teff = l E " , Leff~ 165 km. (9) 
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L e ff corresponds to the distance between KamLAND and the closest set of reactors which 
provides the large fraction of the antineutrino flux. Notice that L e ff is smaller than the 
average (weighted with power) distance 189 km. 

Let us consider the 95% allowed regions. 

(i) The lowest "island" allowed by KamLAND with Am 2 < 2 • 1CT 5 eV 2 corresponds to 
the oscillation phase e // < rc/2. This region is excluded by the absence of significant day- 
night asymmetry of the Super-Kamiokande signal. In this domain, the predicted asymmetry 
at SNO, A S D N N ° > 17%, is still consistent with data. 

(ii) The second allowed region, Am 2 = (5 — 10) • 1CT 5 eV 2 , corresponds to the first 
oscillation maximum, <p e ff ~ ir (maximum of the survival probability). It contains the best 
fit point. 

(iii) The third island is at Am 2 = (13 — 23) • 10~ 5 eV 2 : it corresponds to the oscillation 
maximum (second maximum of the survival probability) with (f) e ff ~ 2n. 

There is a continuum of the allowed regions above Am 2 ~ 3 • 1CT 4 eV 2 . The third region 
merges with the continuum at 99% CL. 

At the la level the second island is the only allowed region. 

2) KamLAND rate. In fig.l we show the regions excluded by the KamLAND rate at 95% 
C.L.. The borders of these regions coincide with contours of constant R 7 ^ = 0.80 and 
B™£ = 0.42 obtained in [16]). The exclusion region at Am 2 - (2 - 5) • 10" 5 eV 2 corre- 
sponds to the first oscillation minimum (minimum of the survival probability) at KamLAND 
(<f) e ff ~ 7r/2). Here the suppression of the signal is too strong. Another region of significant 
suppression (second oscillation minimum with <p e ff = 3n/2) is at Am 2 ~ (9 — 12) • 10~ 5 eV 2 . 

3) . KamLAND spectrum. Low threshold. We analyze the spectrum with the threshold 
Ep h = 0.9 MeV, using the same prescription for the contribution of the geological neutrinos 
as in [1]. The best fit point shifts with respect to (8) to larger mixing and smaller mass 
squared difference: 

Am 2 = 6.86 • 10~ 5 eV 2 , tan 2 9 = 0.48 . (10) 

In fig. 2 we show the allowed regions at 68 %, 95% and 99% confidence levels. 

With lowering the threshold the sensitivity to the spectrum shape increases: the mea- 
surements become sensitive not only to the dominant oscillation peak but also to the lower 
energy oscillation maximum. As a result, the analysis leads to stronger exclusion of the 
oscillation parameter space. In particular, maximal and large mixing parts become less 
favored than in the case of high threshold. 
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4). KamLAND spectrum and rate. Following procedure in [1] we have performed also 
combined analysis of spectrum and rate introducing the free normalization parameter of the 
spectrum, Rkl, and defining the x 2 as 




(11) 



where 




Rkl - 0-611 
0.094 



2 



(12) 



We find results which are very close to those from our spectrum analysis. In particular, the 
best fit value of mixing is tan 2 # = 0.48 and Am 2 = 7.31 • 10~ 5 eV 2 . 

As it follows from our consideration here, the values of oscillation parameters extracted 
from the KamLAND data, (Am 2 , tan 2 9) KL) are in a very good agreement with the values 
from independent solar neutrino analysis (Am 2 , tan 2 #) snn , For the best fit points (1), (8) 
(10) we conclude that within la (see fig. 1 and fig. 2) 



At the same time, the data do not exclude that the solar and KamLAND parameters are dif- 
ferent, and moreover, the difference still can be large. For instance, (Am 2 , tan 2 Q)kl can co- 
incide with the present best fit point or be in the high Am 2 island, whereas (Am 2 , tan 2 6) sun 
can be at lower Am 2 . 

3 Solar Neutrinos 

We use the same data set and the same procedure of analysis as in our previous publica- 
tion [16]. Here the main ingredients of the analysis are summarized. 

The data sample consists of 

- 3 total rates: (i) the Ar-pro duct ion rate, Qati from Homestake [5], (ii) the Ge— production 
rate, Q Ge from SAGE [7] and (iii) the combined Ge— production rate from GALLEX and 



- 44 data points from the zenith-spectra measured by Super-Kamiokande during 1496 days 
of operation [9] ; 

- 38 day-night spectral points from SNO [13, 14]. 

Altogether the solar neutrino experiments provide us with 81 data points. 

All the solar neutrino fluxes, but the boron neutrino flux, are taken according to SSM 
BP2000 [17]. The boron neutrino flux is treated as a free parameter. For the hep— neutrino 
flux we take fixed value F hep = 9.3 x 10 3 cm' 2 s" 1 [17, 26] . 



(Am 2 , tan 2 #)^ L = (Am 2 , tan 2 #) 



sun- 



(13) 



GNO [8]; 
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Thus, in our analysis of the solar neutrino data as well as in the combined analysis of 
the solar and KamLAND results we have three fit parameters: Am 2 , tan 2 # and 

We define the contribution of the solar neutrino data to x 2 as 

Xlun = xlate + XSK + XsNO' (14) 

where Xrata X%k an d xIsno are the contributions from the total rates, the Super-Kami okande 
zenith spectra and the SNO day and night spectra correspondingly. The main result of 
analysis performed in [16] is given here in Eq. (1). 

4 Solar neutrinos and KamLAND 

We have performed two different combined fits of the data from the solar neutrino experi- 
ments and KamLAND. 

1) KamLAND rate and solar neutrino data. There are 81 (solar) + 1 (KamLAND) data 
points - 3 free parameters = 79 d.o.f.. We define the global x 2 for this case as 

X 2 S un,R = Xlun + Xri (15) 

where Xsun ano - Xr are given in (14) and (12). The minimum x^ ni? (min)/(i.o./. = 65.2/79 
corresponds to the C.L. = 86.7% . It appears at 

Am 2 = 6.03 • 10" 5 eV 2 , tan 2 9 = 0.411, f B = 1.05. (16) 

This point practically coincides with what we have obtained from the solar neutrino analysis 
only. 

We construct the contours of constant confidence level in the (Am 2 — tan 2 9) plot using 
the following procedure. We perform minimization of xlun r with respect to for each point 
of the oscillation plane, thus getting x\un r^™ 2 ■> t an2 Then the contours are defined by 
the condition Xstm,fl(Ara 2 , tan 2 9) = xLn,i?(™ n ) + ^X 2 -, where A% 2 = 2.3, 4.61, 5.99, 9.21 
and 11.83 are taken for la, 90%, 95% and 99% C.L. and 3o\ The results are shown in fig. 3. 

According to the figure, the main impact of the KamLAND rate is strengthening of 
bound on the allowed region from below due to strong suppression of the KamLAND rate 
at Am 2 = (2 — 5) • 10~ 5 eV 2 (see fig. 1) - region of the first oscillation minimum in Kam- 
LAND. The lines of constant confidence level are shifted to larger Am 2 . The KamLAND 
rate leads to a distortion (shift to smaller mixings) of contours at Am 2 ~ 10~ 4 eV 2 where 
the second oscillation minimum at KamLAND is situated. The upper part of the allowed 
region is modified rather weakly. 
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2). KamLAND spectrum and solar neutrino data. We calculate 

Xglobal Xsun Xspec 5 

(17) 

where x%ec nas been defined in (7). In this case we have 81 (solar) + 13 (KamLAND) data 
points - 3 free parameters = 91 d.o.f.. The absolute minimum, X 2 giobaA m ^ n ) = 68.2 (which 
corresponds to a very high confidence level: 96.3%), is at 

Am 2 = 7.32 • l(T 5 eV 2 , tan 2 9 = 0.409, f B = 1.05. (18) 

The best fit value of Am 2 is slightly higher than that from the solar data analysis. The solar 
neutrino data have higher sensitivity to mixing, whereas the KamLAND is more sensitive 
to Am 2 , as a result, in (18) the value of Am 2 is close to the one determined from the 
KamLAND data only (8), whereas tan 2 9 coincides with mixing determined from the solar 
neutrino results (1). 

We construct the contours of constant confidence levels in the oscillation plane, similarly 
to what we did for the fit of the solar data and the KamLAND rate (fig. 4). 

As compared with the solar data analysis, KamLAND practically has not changed the 
upper bound on mixing, but strengthened the bound on Am 2 . At the 3<r level we get: 

Am 2 < 2.8 • 10" 4 eV 2 , tan 2 9 < 0.84, 99.73% C.L. . (19) 



The spectral data disintegrate the LMA region. At the 3<r level only a small spot is left 
in the range Am 2 > 2.5 • 10~ 4 eV 2 . 

At the 99% C.L. the rest of the region splits into two "islands" which we will refer to 
as the lower (1-) and higher (h-) LMA regions. (Existence of these two islands can be seen 
already from an overlap of the solar and KamLAND allowed regions in [1]). The 1-region is 
characterized by 

Am 2 = (4.7- 10) • 10~ 5 eV 2 , tan 2 9 = 0.27 - 0.75. (20) 
It contains the best fit point (18). The h- region is determined by 

Am 2 = (12 - 20) • 10~ 5 eV 2 , tan 2 9 = 0.29 - 0.63, (21) 
with best fit point at 

Am 2 = 14.5 • 10~ 5 eV 2 , tan 2 9 = 0.41. (22) 

This point is accepted with respect to the global minimum (18) at about 2a level. 

The excluded region between two islands corresponds to the second oscillation minimum 
(4> e ff ~ 37r/2) at E p ~ (3 — 4)MeV which contradicts the spectral data. 
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Features of spectrum distortion. In the fig. 5 we show the prompt energy spectra of events 
for the best fit points in the 1- and h-regions, Ni(E p ) and Nh(E p ). The spectra can be well 
understood in terms of the effective oscillation phase 4> e ff- 

N(E p ) ~ N (E p ) [l + D{E P ) sin 2 e// ] , (23) 

where D(E P ) is the averaging factor. 

In the best fit point of 1-region (18), the peak at E p « 3.6 MeV corresponds to the 
oscillation maximum <p e ff = it (exact position of maximum of the survival probability is at 
E p = 4.3MeV). The closest oscillation minimum (phase e // = 37r/2) is at E p m 2.4 MeV 
and the next maximum (<f> e ff = 27r) is at £ p « 1.8 MeV. Due to strong averaging effect the 
structures below (in E p ) the main maximum are not profound and look more like a shoulder 
below the peak. The first oscillation minimum is at E p = 7.2 MeV. 

In the case of h-region (22) the spectrum has similar structure but (due to larger Am 2 ) 
the energy intervals between maxima and minima decrease. The main peak at E p m 3.4 
MeV corresponds now to the second oscillation maximum with <p e ff = 2n. One can calculate 
then that the next minimum ( (f) e ff = 5ir/2) is at E p R3 2.7 MeV and the next maximum 
((j> e ff = 37r) is at E — 2.2 MeV. The higher energy oscillation minimum (0 e // = 37r/2) is at 
E = 5 MeV. This pattern can be seen in the fig. 5 which confirms validity of the effective 
phase consideration. 

The measured spectrum, indeed, gives a hint of existence of the low energy shoulder. 
Evidently with the present data it is impossible to disentangle the 1- and h- spectra. Substan- 
tial decrease of errors is needed. Also decrease of the energy threshold will help. According 
to fig. 5, Ni(Ep) > N h (E p ) at E p > 3.5 MeV, and Ni(E p ) < N h (E p ) at lower energies, espe- 
cially in the interval E p = (2.0 — 2.5) MeV. Therefore for the low threshold the difference 
between Ni(E p ) and Nh(E p ) can not be eliminated by normalization (mixing angle). See 
similar discussion in [19]. 



Pull-off diagrams. To check the quality of the fit we have constructed the pull-off dia- 
grams [32] which show deviations, Dk, of the predicted values of the observables K^f from 
the central experimental values, K exp , in the units of the la experimental errors, o~k,' 

D K = K bf- K ^P ^ (24) 
a K 

where K = Q Ar , Q Ge , CC/NC, A S D N N °, A S D K N , R KL , f B . 

In fig. 6 we show the pull-off diagrams for the best fit points from the 1- and h-regions. In 
both regions the largest deviation is for the Ar-production rate. In the h-region: Q Ar > 3.1 
SNU, which is 2.5(7 above the Homestake result. In the 1-region the Qa t deviation is smaller: 
about 1.8(7. 

Other deviations are about la or smaller; they are even smaller in the 1-region. In the 
h-region a very small D-N asymmetry is expected. The 1- and h-regions lead to the opposite 
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sign deviations for the CC/NC ratio at SNO, and therefore future precise measurements of 
this ratio will discriminate among 1- and h- solutions. Notice also that in h-region f B < 1. 

5 Effect of 1-3 mixing 

We assume that Am 2 13 = Am 2 atm = (2 - 3) • KT 3 eV 2 , and there is a non-zero admixture 
of v e in the third mass eigenstate described by the angle 9\ 3 . Both for KamLAND and for 
solar neutrinos the oscillations driven by Am 2 3 are averaged out and signals are determined 
by the survival probabilities 

P ee = (1 - sin 2 9 13 ) 2 P 2 + sin 4 13 « (1 - 2 sin 2 8 13 )P 2 . (25) 

Here P 2 = P 2 (Am 2 2 , #12) is the two neutrino vacuum oscillation probability for KamLAND 
and it is the two neutrino conversion probability for solar neutrinos. The factor (1 — sin 2 9i 3 ) 2 
leads to additional suppression of the KamLAND signal and shift of the allowed regions to 
smaller 9 12 (see detailed study in [27]). The effect of 9 13 on the solar neutrino analysis has 
been discussed in [16]. 

We have performed the combined analysis of the solar neutrino and KamLAND data 
for a fixed value of sin 2 #13 taken at the upper bound given by the CHOOZ experiment [28]: 
sin 2 #13 = 0.04. The number of degrees of freedom is the same as in the 2u fit and we follow 
procedure described in sect. 4 with survival probabilities modified according to (25). 

In fig. 7 we show results of the combined analysis of the solar data and the KamLAND 
rate. As can be concluded from comparison of fig. 3 and fig. 7, the main impact of non- 
zero 1-3 mixing is a slight shift of the allowed region to larger Am 2 and smaller tan 2 9. In 
particular, we find that the best fit, Xmin/d-O-f- = 66.21/79, is at 

Am 2 = 6.60 • 10" 5 eV 2 , tan 2 9 = 0.411 f B = 1.08. (26) 

In fig. 8 we show result of the combined analysis of the solar neutrino data and the 
KamLAND spectrum. From fig. 8 and fig. 4 we observe the same effect: 1-3 mixing leads 
to a shift of the allowed region to larger Am 2 and smaller tan 2 9. In the best fit point we 
find xLJd-o.f. = 69.24/91 and 

Am 2 = 7.17- 10~ 5 eV 2 , tan 2 9 = 0.41, f B = 1.07. (27) 

The introduction of the 1-3 mixing slightly worsen the fit: A% 2 ps 1, it requires higher 
original boron neutrino flux, than in the 2v- case. At the same time, the 1-3 mixing improves 
a fit in the h-region: It is accepted now at 90 % C.L. with respect to the best fit point (27). 
This region merges with the small spot at high Am 2 . The best fit point in the h-region is 
shifted to smaller mixing: tan 2 9 = 0.36. 

For convenience, the results of different fits are summarized in the Table I. It shows 
high stability of the extracted parameters with respect to a type of analysis. 
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Type of data fit 


A 2 

Amr 


tarn 9 


fB 


9 / 1 f 
Xmm/cl.O.f. 


solar 2v 


6.15 


0.406 


1.05 


65.2/78 


KL spectrum, 2u 


7.34 


0.453 




2.92/11 


solar + KL rate, 2v 


6.03 


0.411 


1.05 


65.3/79 


solar + KL spectrum, 2v 


7.32 


0.409 


1.05 


68.2/91 


solar + KL rate, 3v 


6.60 


0.411 


1.08 


66.2/79 


solar + KL spectrum, 3u 


7.17 


0.409 


1.07 


69.2/91 



Table 1: The parameters of the best fit points as well as Xmin/ 'd.o.j '. from different analyses 
of the data; Am 2 is in the units 10~ 5 eV 2 . 



6 Next step 

The key problems left after the first KamLAND results are 

• more precise determination of the neutrino parameters: in particular, (i) precise de- 
termination of the deviation of 12-mixing from maximal mixing, (ii) strengthening of 
the upper bound on Am 2 (iii) discrimination between the two existing regions; 

• searches for effects beyond the single Am 2 and single mixing approximation; 

• searches for differences of the neutrino oscillation parameters determined from Kam- 
LAND and from the solar neutrino experiments. 

Notice that precise knowledge of the parameters is crucial not only for the neutrinoless 
double beta decay searches, long baseline experiments, studies of the atmospheric and super- 
nova neutrinos, etc., but also for understanding of physics of the solar neutrino conversion. 
In the region of the best fit point a dominating process (at least for E > (0.5 — 1) MeV) 
is the adiabatic neutrino conversion (MSW), whereas in the high Am 2 region allowed by 
KamLAND at the 2a level, the effect is reduced to the averaged vacuum oscillations (a la 
Gribov-Pontecorvo) [29] with small matter corrections. 

In this connection, we will discuss two questions. 

How small Am 2 can be? This is especially important question, e.g., for measurements of 
the earth regeneration effect. In the 1- region we get 

Am 2 > 4 • 10~ 5 eV 2 , (3a) (28) 

and it is difficult to expect that lower values will be allowed. 
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The bound (28) appears as an interplay of both the KamLAND rate and the shape. 
As it follows from the fig. 3, the KamLAND rate strengthens the lower bounds: Am 2 > 
(3.0, 4.0, 4.7) x 10~ 5 eV 2 , at the la, 2a, 3a correspondingly which should be compared 
with Am 2 > (2.5, 3.2, 4.0) x 10~ 5 eV 2 from the solar analysis only. Adding the spectral 
data results in the bounds Am 2 > 3.3, 5.5, 6.2 x 10" 5 eV 2 

With decrease of Am 2 the oscillatory pattern of the spectrum shifts to lower energies. 
For Am 2 = 5 • 10~ 5 eV 2 the maximum of spectrum is at E p = 2.7 MeV and the oscillation 
suppression increases with energy [16]. The oscillation minimum is at E p 5 MeV. If 
Rkl(2.7 MeV) = 0.81, then ifo L (4.0 MeV) = 0.47. The KamLAND spectrum does not 
show such a fast decrease. 

One can characterize the spectrum distortion by a relative suppression of signal at the 
high (say, above 4.3 MeV) and at the low (below 4.3 MeV) energies 1 . The energy interval 
(2.6 - 4.3) MeV contains KamLAND energy 4 bins. Introducing the suppression factors 
Rkl{< 4.3 MeV) and Rrl{> 4.3 MeV) we can define the ratio 

= 1 — Rkl(> 4-3 MeV) 

1 - R K l{< 4.3 MeV) V ; 

k which we will call the shape parameter does not depend on the mixing angle and nor- 
malization of spectrum. It increases with increase of the oscillation suppression at high 
energies. 

Using the KamLAND data we get the experimental value 

k exp = 0.84+ 42 5 , la. (30) 



In fig. 9 we present the dependence of the shape parameter k on Am 2 for fixed mixing: 
tan 2 9 = 0.41. For the spectrum which corresponds to the best combined fit we find k = 0.70, 
whereas for Am 2 = 5 ■ 10~ 5 eV 2 the ratio equals k = 2.0. For the h-region best point: 
k = 0.94. 

Notice that in the 1-region below Am 2 = 8 • 10~ 5 eV 2 , k increases quickly with decrease 
of Am 2 , reaching a maximal value at Am 2 = 5.5 ■ 10~ 5 eV 2 . Below that, the parameter k 
decreases with Am 2 . In this region, however, the total event rate decreases fast giving the 
bound on Am 2 . This explains a shift of the allowed (at 3a) region to smaller mixings with 
decrease of Am 2 . 

Notice also that the central experimental value of k can be reproduced in the both al- 
lowed regions (1- and h-). Therefore future precise measurement of spectrum will further 
sharpen determination of Am 2 within a given island. To discriminate among the islands 
one needs to use more elaborated criteria (not just k) or a complete spectral information. 



An alternative way is to introduce the moments of the energy spectrum [30] . 
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How large is the large mixing? In contrast to [1] our best fit point is at non- maximal mixing 
(sin 2 29 = 0.86) being rather close to the best fit point from the solar neutrino analysis. 
Similar deviation from maximal mixing has been obtained in our rate+spectrum analysis. 
Notice that the KamLAND data have weak sensitivity to the mixing (weaker than the solar 
neutrino data). The allowed regions cover the interval 

tan 2 9 = 0.12 - 1.00 (0 < tt/2), 95%C.L.. (31) 

Even at la the interval tan 2 9 = 0.23 — 1.00 (9 < n/2) is allowed. The reason is that 
KamLAND is essentially the vacuum oscillation experiment (see evaluation of matter effects 
in KamLAND in [23]), and effects of the vacuum oscillations depend on deviation from 
maximal mixing which can be characterized by e = (1/2 — sin 2 9) quadratically: P oc 1 — 4e 2 . 
The matter conversion depends on e linearly: P oc 1 — 2e [31]. 

Notice that maximal mixing is rather strongly disfavored by all measured solar neutrino 
rates. In the point indicated by KamLAND we predict the charged current to neutral 
current measurement at SNO, the Argon production rate and the Germanium production 
rate: 

CC 

— = 0.51 (+3.3(j), Q Ar = 3.2 SNU (+2.5a), Q Ge = 63 SNU (-1.6a). (32) 

IN 

In brackets we show the pulls of the predictions from the best fit experimental values. 
As follows from the fig. 10 future precise measurements of the CC/NC ratio at SNO will 
strengthen the upper bounds on both mixing and Am 2 . 



It is important to "overdetermine" the neutrino parameters measuring all possible ob- 
servables. This will allow us to make cross-checks of selected solution and to search for 
inconsistencies which will require extensions of the theoretical context. In this connection 
let us consider predictions for the forthcoming measurements. 

1). Precise measurements of the CC/NC ratio at SNO. In fig. 10 we show the contours of 
constant CC/NC ratio. We find predictions for the best fit point and the 3a interval: 

CC 

— = 0.33^, 3a. (33) 

At 99% C.L. the islands split and we get the following predictions for them separately: 
CC CC 

^(1 - region) = 0.33^ , ^(h - region) = 0.38±° : °g , 99% C.L. . (34) 

Values of CC/NC < 0.35 will exclude the h-region. Precise measurements of the ratio 
CC/NC will also strengthen the upper bounds on mixing and Am 2 . 
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2). The day-night asymmetry at SNO. The KamLAND provides a strong lover bound on 
Am 2 , and shifts the best fit point to larger Am 2 . This further diminishes the expected 
value of day- night asymmetry. In fig. 10 we show the contours of constant Ap%°. The best 
fit point prediction and the 3a bound equal 

A S N N D ° = 2.8±0.8% , (la), A S N N ° < 9% (3a). (35) 

The present best fit value of the SNO asymmetry, 7%, is accepted at about 99% C.L.. Obser- 
vations of the asymmetry Aff^f > 1% will exclude the h-region. The expected asymmetry 
at Super-Kamiokande is even smaller: In the best fit point we expect Aff£ (1.7 — 2.0)%. 



3). The turn up of the spectrum at SNO and Super-Kamiokande at low energies. Using 
results of [32] we predict for the best fit point (18) an increase of ratio of the observed to 
expected (without oscillation) number of the CC events, Rsnoi from 0.31 at 8 MeV to 0.345 
at 5 MeV, so that 

In Super-Kamiokande the turn up is about (5 — 7)% in the same interval (5 - 8) MeV. 



4). Further KamLAND measurements. Possible impact can be estimated using figures 1, 
2, 5. A decrease of the error by factor of 2 (which will require both significant increase of 
statistics and decrease of the systematic error) will allow KamLAND alone to exclude all the 
regions but the 1-region at 95 % C.L., if the best fit is at the same point as it is determined 
now. 



5). BOREXINO. At 3<x we predict the following suppressions of signals with respect to the 
SSM predictions to BOREXINO experiment [33], in the two allowed regions: 

R B (\ - region) = 0.61 -0.73, R B (h - region) = 0.62 - 0.73. (37) 

So, BOREXINO will perform consistency check but it will not distinguish the 1- and h- 
regions. 



6). Gallium production rate. In the best fit point one predicts the germanium production 
rate Qc e = 71 SNU. In the h-region best fit point Qc e = 72 SNU. So, the Gallium experi- 
ments do not discriminate among the 1 - and h - regions. However, precise measurements of 
Qce are important for improvements of the bound on the 1-2 mixing and its deviation from 
maximal value. 



13 



7 Conclusions 



1. The first KamLAND results (rate and spectrum) are in a very good agreement with the 
predictions based on the LMA MSW solution of the solar neutrino problem. 

2. Our analysis of the KamLAND data reproduces well the results of the collaboration. 
The oscillation parameters extracted from the KamLAND data and from the solar neutrino 
data agree within la. 

3. We have performed a combined analysis of the solar and KamLAND results. The main 
impact of the KamLAND results on the LMA solution can be summarized in the following 
way. KamLAND 

• shifts the Am 2 to slightly higher values 6.15 • 1(T 5 -> 7.3 • 10" 5 eV 2 ; (the mixing is 
practically unchanged: tan 2 9 = 0.41, and this number is rather stable with respect to 
variations of the analysis; 

• establishes rather solid lower bound on the Am 2 : Am 2 > 4 • 10~ 5 eV 2 (3a); 

• disintegrates the LMA region at 99% C.L. into two parts. 

KamLAND further disfavors high values of Am 2 : Am 2 > 3 • 10~ 4 eV 2 . 

4. Inclusion of the 1-3 mixing shifts the allowed regions to larger Am 2 and smaller tan 2 9. 
X 2 slightly increases with sin 2 #i 3 . 

5. The KamLAND results strengthen the upper bound on the expected value of the day- 
night asymmetry at SNO: A s ^f < 9%. We predict about 10% turn up of the energy 
spectrum at SNO in the interval of energies (5 - 8) MeV. The CC/NC ratio is expected to 
be CC/NC ~ 0.33, that is, near the present best SNO value. 

Future SNO measurements of the CC/NC ratio and Aj^ will have further strong 
impact on the LMA parameter space. 
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Figure 1: The KamLAND spectrum analysis for E p > 2.6 MeV. Shown are the allowed 
regions of oscillation parameters at 68% (inner solid lines), 95% (grey) and 99% C.L. (outer 
solid lines). The best fit point is indicated by star. Also shown are the regions excluded by 
the rate analysis at 95% C.L. (dark). 
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Figure 3: The allowed regions in tan 2 9 — Am 2 plane from a combined analysis of the solar 
neutrino data and the KamLAND rate, at la, 90%, 95%, 99% and 3<r C.L.. The best fit 
point is marked by star. 
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Figure 4: The allowed regions in tan 2 9 — Am 2 plane, from a combined analysis of the solar 
neutrino data and the KamLAND spectrum at la, 90%, 95%, 99% and 3cr C.L.. The best 
fit point is marked by star. 
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Figure 5: The expected prompt energy spectra for the best fit points from the 1- region (solid 
histogram) and h- region (dashed histogram). Also shown are the KamLAND experimental 
points. 



21 



NC/CC 




Figure 6: The pull-off diagrams for the best fit points in the 1- and h- regions. Shown 
are deviations of the predicted values of different observables from the central experimental 
values in the units of la (experimental). 
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Figure 7: Three neutrino analysis with sin 2 #13 = 0.04. The allowed regions in tan 2 6 — Am 2 
from a combined fit of the solar neutrino data and the KamLAND rate at the la, 90%, 
95%, 99% and 3a C.L.. 
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Figure 8: Three neutrino analysis with sin 2 #13 = 0.04. The allowed regions in tan 2 6 — Am 2 
plane from a combined analysis of the solar neutrino data and the KamLAND spectrum at 
the la, 90%, 95%, 99% and 3a C.L.. 
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Figure 9: The dependence of the shape parameter k on Am 2 for tan 2 6 = 0.41. Shown are 
the central experimental value (dotted line) and the la experimental band (shadowed). 
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Figure 10: Predictions for the CC/NC ratio and the Day-Night asymmetry at SNO. The 
dashed lines are the lines of constant CC/NC ratio (numbers at the curves) and the dotted 
lines show the lines of constant A'^f (numbers at the curves in %). We show also the 
(la and 3<r) allowed regions of the oscillation parameters from the combined fit of the solar 
neutrino data and the KamLAND spectrum. The best fit point is indicated by a star. 
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